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Abstract. We consider m two-dimensional semi-infinite planes of Ising spins joined 
together through surface spins and study the critical behaviour near to the junction. 
The m = limit of the model — according to the replica trick — corresponds to the 
semi-infinite Ising model in the presence of a random surface field (RSFI). Using 
conformal mapping, second-order perturbation expansion around the weak— and 
strong-coupled planes limits and differential renormalization group, we show that the 
surface critical behaviour of the RSFI model is described by Ising critical exponents 
with logarithmic corrections to scaling, while at multiple junctions (m > 2) the 
transition is of first-order. There is a spontaneous junction magnetization at the 
bulk critical point. 



The critical behaviour of systems near a plane where translational invariance is 
broken, is of considerable recent interest [1, 2]. The prototype of these problems 
is represented by the critical phenomena at a free regular surface (semi-infinite 
criticality) , and in more complex problems the effect of a perturbation (e.g. surface 
coupling enhancement, interfaces, defects, random surface fields etc.) can be analysed 
by relevance-irrelevance type criteria [3-5]. Such a stability analysis, however, does 
not work for the two-dimensional Ising model in the case of marginal perturbations 
caused by a defect line [3, 4] or a random surface field (rSFI) [5], when the defect 
exponent y^ = 0. In the former case non-universal critical behaviour was found by 
an exact calculation [6, 7], while for the RSFI model no definite answer is known yet. 
One of our purpose in the present Letter is to clarify the critical behaviour of the RSFI 
model. 

To study this problem we introduce a series of models consisting of m semi-infinite 
planes of Ising spins where the spins at different surfaces are joined together by nearest 
neighbour couplings (see figure 1(a)). (We note that in a recent paper Indekeu and 
Nikas [8] introduced a junction as a product of surface spins and studied the wetting 
phenomena in the m = 3 system in the frame of Landau theory) . The Hamiltonian of 
the system is given by: 
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a) b) 

Figure 1. (a) m = 3 semi-infinite planes of Ising spins joined together by surface 
spins and (6) the corresponding system in the strip geometry. 



where 

oo oo 
y= — oo x=l 

and 

oo 

-I3V = J (3) 

y=-oop<p' 

Here a"P(j/, a;) = ±1 are Ising spins at position j/, x on the p-th. plane, and (3 = 1/k^T. 

For m = 1 and m = 2 we obtain a semi-infinite system and two semi-infinite 
systems joined by a defect line, respectively, while for m > 3 as a possible physical 
realization one can imagine interacting magnetic ions segregated along various planar 
grain boundaries, three of which meet along a linear junction. Due to pair interactions 
in the junction (equation (3)) the perturbation represented by V is marginal in the 
ordinary surface transition point for al\ m ^ 1\ thus to clarify the actual critical 
behaviour one needs detailed investigations. 

Now we show that the multiple junction problem is connected to the RSFI model 
defined by the Hamiltonian: 

oo 
y=-oo 

where the random surface field h{y) has a Gaussian distribution: 

Since the disorder is quenched, it is the free energy rather than the partition function 
which must be averaged. Using the replica trick: (log Z) = lim,^^^ [((Z"^) — l)/m] one 
can easily show that the effective Hamiltonian of the problem is just the m— >0 limit 
of equations (l)-(3) with J = A^. 



(4) 



(5) 
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In general we are interested in the critical behaviour near to the junction, e.g. 
we look for the critical exponent rj^ describing the decay of spin-spin correlations 
((7^(0, 1)) oc when the system is at the bulk critical point. Now 

wc suppose that critical correlations in the model transform covariantly under a 
conformal transformation. It is exactly known for m — 1 and m — 2 [9], on the 
other hand the gap-exponent relation (6) might be valid, even if the system is 
not conformally invariant [10]. In the following we map the system onto the strip 
geometry [11], where the calculation is usually simpler to perform. Denoting the 
points on the pth plane by the complex number z^, then the conformal transformation 
Wp = (L/7r) log;2p, p = l,2,...,m, maps the semi-infinite planes onto strips of width L, 
and the surface spins at both ends of the strips are connected to each other with the 
same type of coupling as in the plane geometry equation (3) (see figure 1(b)). 

The critical exponents in the strip geometry can be calculated from the finite-size 
behaviour of the correlation length [11]. More precisely we consider the extreme- 
anisotropic limit [12] of the model, when the transfer matrix along the strip is expressed 
as T = exp(— a^f), where a is the lattice spacing and if is a quantum Hamiltonian. 
Then, following Cardy's derivation [13], one can show that the spectrum of the critical 
Hamiltonian operator H in the large-L limit is given as: 

E.-E^ = ^v^x, (6) 

where Eq and E^ are the ground state and the ith excited state of H, respectively, and 
V,. is a normalizing factor, the so-called sound velocity. The set of critical dimensions 
X- describes the decay of correlations of scaling operators along the junction in the 
plane geometry: {^\0)^^{z)) o<:\z\~^^\ For the spin operator we have 77^ = 20;^. The 
spectrum in (6) usually has a tower-like structure; the levels in the same tower differ 
by an integer from the lowest one: x^^ = + 1; Z = 0, 1, 2, . . . , and x^ is the critical 
dimension of a primary operator [11]. 

For our model, the quantum Hamiltonian is given as 

m 

L-l L 

V=-^iYl -^lYI <^L)ai{L) (7b) 

p<p' p<p' 

In (7a, h) and are Pauli matrices at chain p on site x, the bulk critical point 
corresponds to h = 1 [12], while in the units used in (7) f g = 2 [14]. Note that we 
put different values of the couplings at both ends of the chains (Ai,A^): the model in 
(1-3) corresponds to X^ = Xj^ = X. 

We mention that one may also consider systems composed of planes having two free 
surfaces, say at a; = and y — 0. For these ^half-infinite^ systems, in the Hamiltonian 
equations (2), (3), the summations over y run from to 00. Now in the transformed 
geometry the strips are coupled only at one edge and the others are free; furthermore 
in the transformation tt is replaced by 7r/2, the angle at the corner [11]. In this case 
in (7b) we have Aj^ = A and A^ = 0. 
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To calculate the spectrum of (7) one has to keep in mind the possible presence 
of strong logarithmic corrections to scaling, and therefore to try to push analytical 

calculations as far as possible. Since for general value of m no exact solution can 
be found, we perform a perturbation expansion around the uncoupled chains limit 
and sum up the most diverging contributions by using the differential renormalization 
group technique. 

The actual calculation of the coeflBcients of the weak-junction expansion is rather 
cumbersome; therefore here we present only the final result of the second-order 
calculation; details of the derivation will be published elsewhere [15]. The first gap of 
the system — corresponding to magnetic excitations — in the large-L limit is given by: 



\-{m- + A^) + (m - l)^AiA^ 



_ 1)(^ _ 2) A log L-Cj {Xj + Ai) + 



+ 0{1/L') (8) 



where C = 0.03672... is a constant, and according to (6) the quantity in the square 
bracket is just the leading magnetic exponent x^. The main observation concerning 
equation (8) is that in first order is regular and coupling dependent, but the second- 
order coefficient diverges as log L. The perturbation series for higher gaps shows the 
same qualitative picture. On the basis of the differential renormalization group [16] 
one assumes that the higher-order terms of the expansion in (8) are divergent too, 
but these singular terms sum up to a regular contribution when the perturbation is 
marginally irrelevant. 

To show this we start writing the differential renormalization group equations of 
the problem under a change of the length scale e' in the form [16]: 

^=y^X,+bXl + 0{Xl) 

(9) 

^ = y,^L + bxl + o(Ai) 

where the defect exponent y^ = Q and the perturbation is marginally irrelevant and 
marginally relevant for b<0 and 6>0, respectively. The solutions of (9) are given as 
X-^^{1) = X^/ (l — bX^l) and A^(Z) = A^/(l — bA^^Z), respectively. From the transformation 
form of the inverse correlation length 

C'(Ai, A^, L-i) = e-^C'(Ai(0, A^(/), L'^eO 
one obtains the scaling prediction of the gap for finite systems: 

Expanding $^ up to second order in powers of A^^, A^ one can verify that its form is 
compatible with the expansion in (8) with 6 = (2/7r)(m — 2). One can verify similarly 
that the scaling form (10) is also valid for higher gaps with the same value of b. This 
fact is due to the structure of the second-order degenerate perturbation calculation: 
different gaps are represented by different secular matrices, but all the matrix elements 
are the same function of A^^, A^ and L, independently of the matrix. 
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Thus we arrive at the conclusion, that the relevance-irrelevance behaviour of the 
perturbation caused by the junction depends on the value of m: for positive defect 
couplings it is marginally irrelevant and marginally relevant for m < 2 and m > 2, 
respectively. In the borderline case m = 2 the perturbation is fully marginal, the 
critical exponents are coupling dependent [6, 7]. 

For m < 2 in the large-L limit: l^jA^logL ^ 1, |6|A£^logL ^ 1, the finite-size 
behaviour of the magnetic exponent is given by: 

1 m - 1 1 

"^=2-23^1^ + - ^''^ 

Concerning the surface transition of the RSFI model (m = 0) we obtain Ising critical 
exponents with logarithmic correction to scaling 



which is one of the main results of our paper. The finite-size form of the exponent Xg 
in (11a) gives a theoretical basis to deduce effective exponents from MC simulations or 
from transfer matrix calculations. We note, that similar behaviour — exponents of the 
pure model with logarithmic corrections — has been observed for the magnetization 
and the susceptibility of the random-bond two-dimensional Ising model [17], which 
can be considered as the bulk analogue of the RSFI model. 

Now we turn to the problem of multiple junctions (m > 2) with ferromagnetic 
interactions. In this case the perturbation is marginally relevant, consequently the 
critical behaviour is controlled by a new fixed point. We believe that the system for 
non-zero defect couplings undergoes a first-order transition, there is a spontaneous 
junction magnetization at the bulk critical point. To prove this conjecture we consider 
the strong defect limit of the problem and investigate the stability of the corresponding 
fixed point. 

In the limit Aj^ — > oo and/or Aj^ — > oo all the spins in the junctions are parallel, 

consequently the ground state of the system is two-fold degenerate, i.e., the lowest gap 
is zero and there is a spontaneous junction magnetization of the system. In this limit 
the spectrum of H can be constructed as the direct product of m Virasoro algebras 
corresponding to the spectrum of the Ising model at the extraordinary transition 
point [11]. 

The stability of this fixed point is determined through the size dependence of the 
lowest gap for finite defect couplings. In the half-infinite case, Aj^^-l, A£^ = in first 
order of 1/ Aj^ the gap is given by [14] 

1 / 1 A "^^^ 

^1"^°= (m - 7)2-2 Ar^ [l) +••• (12) 

while for non-zero A^ the leading finite-size dependence of the gap remains the same 
with a more complicated prefactor than in (12). 

According to (12) for m > 2 at finite couplings the first gap vanishes faster than 
1/L, thus from (6) Xg = and the perturbation is irrelevant. The ground state of the 
system is degenerate and there is a spontaneous junction magnetization at the bulk 
critical point in accordance with our claim. We note that similar mechanism has been 
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observed for other problems where surface or interface ordering take place at the bulk 
critical point: the lowest gap vanishes algebraically with the size of the system [18-20]. 
For m = 2 according to (12) the perturbation is marginal, is coupling dependent, 
while for m<2 the gap vanishes slower than 1/L, thus the extraordinary transition is 
unstable in this region. 

Finally we turn to discuss the critical behaviour near to a junction with antiferro- 
magnetic (Af) interactions. In this case according to equations (9) the signs of b 
for marginally irrelevant and marginally relevant perturbations are interchanged with 
respect to the ferromagnetic junction. As a consequence for m>2 the perturbation at 
the ordinary surface transition point is marginally irrelevant; the critical exponent 
together with the logarithmic finite size correction is given by (11). 

It is interesting to study the strong junction limit for AF couplings, which gives 
qualitatively different results for m=odd and m=even. For m=odd, due to frustra- 
tion, there is no extraordinary transition even for infinitely strong couplings in the 
defect: the first gap is proportional to L~^/^ [15]. This phenomenon resembles the 
absence of phase transition at T = for super-frustrated models [21]. On the other 
hand for m=even there is an extraordinary transition for infinitely strong AF couplings. 
The perturbation to the first gap is marginal in leading order: Ei—Eq(x1/{XL), thus 
further analysis is needed to decide on the type of transition in this case. 
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